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INTRODUCTION 


Tho aosczlo ^tftorbo'ly is i*>ne of the m^^in drag-produclnij components of ixn aircraft 
propulsion system. Thus, considerable effort has boon devoted to developing toch- 
niquQS for predicting tho afterbody flow field and drag. Using methods ranging from 
empirical techniques based on experimental data (ref. 1) to solutions of the Navier^ 
Stokes equations with turbulence closure (refs. 2 and 3), investigators have had some 
success for isolated axisymmetric nozzles. Patched viscous'-inviscid interaction 
methods (refs. 4 to 6, for example) have had considerable success* The patched 
methods rt'quire considerably less computational time than the Navier-Stokes solutions 
and generally yield results of comparable (and in some cases greater) accuracy. 

Tho present paper describes a viscoua-inviscid method that has evolved over a 
number of years through tho efforts of several researchers. Reubush and Putnam 
(ref. 4) combined a conventional Integral boundary-layer technique (ref. 7) with a 
linearized potential flow computation to account for the boundary- layer dlaplacemoat 
effect. For separated flows, they employed the di scr iminat Ing streamline concept of 
Presz (ref. 8) to separate the reverse flow region from the outer flow. Wilmoth 
(ref. S) extended the method to transonic speeds by replacing the linearized inviscid 
method with the relaxation procedure of South and Jameson (ref* '3) for solving the 
exact nonlinear potential flow equation. Abeyounls (ref. 10) conducted oil-flow 
studies to determine the separation location on a series of circular-arc boattails 
and used the results to evaluate several prediction techniv^ues including the Presz 
model. Presz ot al. (ref. 11) developed an improved analytical model of tho sepa- 
rated region which accounts for axial pressure gradients, surface skin friction, and 
jet entrainment. Putnam (ref. 6) incorporated this improved model together with a 
simple one-dimensional model of the inviscid jet plume hlockvage Into the linearized 
method. Wilmoth and Dash (ref. 12) developed a jet entrainment model that uses the 
overlaid mixing analysis of Dash and Pergament (ref. 13) to calculate a jet plume 
displacement. This entrainment model was incorporated into the transonic method by 
Wilmoth (ref. 14) together with the shock-captur inir— ahock-f itt ing inviscid plume 
model of Dash and Thorpe (ref. 15). Tho work described above has provided the basis 
for the present viscous- inviscid model. 

The present model iteratively combines the South-Jameson relaxation procedure, 
the Resl)otka-Tucker boundary-layer solution, the Presz separation model, the Daah- 
Porgament mixing model, and the Dash-Thorpe inviscid plume model. These computa- 
tional models are combined into a single computer program called RAXJKT* RAXJET 
px'edicts transonic axisymmetric flow over nozzle afterbodies with supersonic jot 
exhausts and accounts for boundary-layer displacement, separation, jet entrainment, 
and inviscid jot plume blockage. This j’japer describes the various components of tho 
computational model, illustrates its capabilities by comparison with experimental 
data, and provides a user's guide to the computer program. The computer program may 
be obtained as program LAR- 12957 from 
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SYMBOLS 


'ref 


The symbols used in the computer printouts are given in the second column, 
cross-sectional area of body of revolution 
SRBP reference area for drag coefficients 

BL boundary layer 

CDP,APT afterbody friction drag coefficient, 

•S(x„) 


•'o/w \ 


C dS 

A J f 

ref •'s(x ) 

A 




CDK,BOD body-of-revolutlon friction dreg coefficient, 

•S(x„) 

dS 


1. f'-H' 


ref *'0 

COP, APT afterbody pressure drag coefficient, 
•A(x„) 


_L_ 

Vef •^A(x^) 
A 


C dA 
P 


COP, BOD body-of-revolution pressure drag coefficient, 
•A(xJ 


CF 


CP 




C dA 
P 


ref •'0 

CDT,APT afterbody total drag coefficient, CDP,AFT CDF, AFT 
CDT,BOD body-of-revolution total drag coefficient, CDP,BOD CDF, BOD 


skin-friction coefficient, t /q 

w » 


pressv^re coefficient, 


P * K 
q 


D 


kt2 


nozzle base diameter 
nozzle maximum diameter 

GAMMA ratio of specific heats of gas mixture 

H boundary-layer shape factor, 

I longitudinal grid index 

two- equation turbulence model 
length of afterbody 
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h 

M 


N 

CD 

NPR 


P 

Pc 

Pt,2 


P 




q 


CO 


r 


S 

T 


u 


V 


X 


L or REPL 

M 

ML 

MO 


P 

PINF 

PT 

PTINP 

R 


RDS 

RKPF 

RPAVG 

RPMAX 

RVI 

T 

TINP 

U 

UINP 

V 

X 


mference length 
Mach numbar 

local surface Mach number 
free-stream Mach n-^mber 

noxsle pressure ratio# ratio of jet total pressure 
to free-stream static pressure 

static pressure 

free-stream static pressure 

total pressure measured by pitot tube; total pressure 
downstream of normal shock if M > 1 

free-stream total pressure 

free-stream dynamic pressure 

radial coordinate of cylindrical coordinate system with 
orlg I at nose of body of revolution 

radius of jet at noesle exit 

body radius corrected for discriminating streamline and 
inviscid plume boundary 

effective body radius used in inviscid external flow calculation 
average residual in relaxation solution of potential equation 
maximum residual in relaxation solution of potential equation 
radius at outer edge of boundary layer or mixing layer 
surface area 
temperature 

free-stream temperature 
axial velocity 
free-stream vela, ity 
radial velocity 

axial coordinate of cylindrical coordinate system with origin at 
nose o. body of revolution 

axial location of start of afterbody 
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^sep 


X,Y 


a 




6 * 

A 

9 


XM dxial length ot body of revolution 

axial location of separation 

computational coordinates of body-fitted grid used in inviscid 
external flow solution 

ALF radial grid stretching parameter used in external flow solution 

boundary- layer or mixing-layer thickness 
DEL* boundary-layer or mixing-layer displacement thickness 

incremental value 

THETA boundary- layer momentum thickness 

body- fitted grid coordinates 
boundary- layer shear stress at body surface 


DESCRIPTION OP METHOD 

The viscous- inviscid interaction model is constructed by dividing the afterbody 
flow field shown in figure 1 into six separate computational regions (see fig* 2)* 

The viscous interaction is taken into account by displacing the body-plume shape with 
an appropriate displacement thickness to yield an effective body geometry for the 
Inviscid external flow calculation* The inviscid external flow solution and inviscid 
jet exhaust solution provide the necessary flow conditions to calculate conditions in 
the viscous regions by parabolized marching procedures* The viscous and inviscid 
flow fields are then iteratively solved* 


Inviscid External Flow Solution 

The inviscid external flow solution is based on the relaxation procedure of 
South and Jameson (ref* 9) for solving the exact nonlinear potential flow equation in 
nonconservative form. This technique Is valid for subsonic, transonic, and super- 
soric Irrotationai flow over bodies of revolution* RAXJET uses the computer program 
RAXBOD developed by Keller and South (ref* 16) to implement this technique. RAXBOD 
usee a body-normal coordinate system from the nose up to the first horizontal tangent 
and a body-fitted sheared coordinate system aft of this point. The computational 
grid is stretched to infinity in both the normal and downstream axial directions to 
faclHtate handling of the far-fleld boundary conditions* Since RAXBOD requires that 
a complete body of revolution be specified, special procedures have been Incorporated 
into ilAXJBT to determine grid stretching parameters that provide optimal grid point 
distribution over the afterbody. 


Inviscid Jet Exhaust Solution 

RAXJET uses the computer program SCIPAC (ref. 17) to solve the flow field In the 
inviscid jet region* SCIPAC solves the flow field by explicit spatial marching of 
the conservative finite-difference form of the inviscid flow equations for a uniform 
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composition gas mixture. The nozzle exit flow is assumed to be supersonic, and the 
calculation is initiated at the exit with exhaust properties prescribed. A uniform 
distribution of grid points is used between the axis and the inviscid plume boundary. 
As the calculation proceeds downstream, oblique shocks are numerically captured while 
the triple point that occurs because of the formation of a Mach disk in an under- 
expanded jet is treated by a shock-f ttting procedure. The Mach disk location is 
determined by the iterative procedure of Abbett (ref. 18) which treats the subsonic 
region downstream of the Mach disk as a one-dimensional isentropic streamtube. 
Generally, for mildly underexpanded jets (low nozzle pressure ratios), the Mach disk 
diameter is small, and the subsonic region is approximated by a constant area 
streamtube. 

For low Mach number (near unity) jets, the supersonic marching procedure becomes 
inefficient since the marching step size of the explicit scheme approaches zero as 
the Mach number approaches one. Also, under such conditions, the shock-fitting pro- 
cedure may yield poor results since shocks are weak and difficult to define from the 
numerical results. Thus, RAX JET includes a low Mach number solution procedure in 
which the flow properties are calculated based on an isentropic expansion to the 
local inviscid external pressure (refs. 15 and 17). This isentropic procedure may 
also be applied to the flow downstream of a Mach disk. The computer program provides 
tests based on exit pressure and Mach number to determine whether the low-Mach-nuraber 
procedure should be applied. 


Boundary-Layer Solution 

The properties in the attached boundary-layer region are solved by a modified 
version of the Pashotko-Tucker integral method for turbulent flows (ref. 7). The 
computer algorithm is described in references 11 and 19. The integral solution 
is obtained by conventional boundary-layer marching procedures to yield the 
displacement-thickness distribution fi*(x) over the body. 

Knowledge of the internal boundary-layer thickness is needed only to initialize 
the jet wake and mixing analysis. Since the nozzle internal boundary layer is gener- 
ally quite thin compared to the afterbody external boundary layer, an exit boundary- 
layer thickness is estimated from flat-plate theory. 


Separated Flow Solution 

The analysis of the separated region consists of predicting the separation loca- 
tion and calculating the discriminating streamline shape. The discriminating stream- 
line analysis is further divided into two regions: the separation region from the 

separation point to the end of the body and the jet wake region from the end of the 
body to the reattachment point (see fig. 2). 

Separation location prediction .- The location of separation on the afterbody is 
calculated by the control volume technique developed by Press (ref* 8)* The tech- 
nique determines the separation point based on a calculated pressure rise to separa- 
tion from the minimum pressure on the afterbody* Abeyounis (ref* 10) con^ared the 
results of the Press model and several other prediction techniques with the results 
of oil-flow studies on a series of circular-arc boattail nossles* While none of the 
methods examined was able to successfully reproduce the experimental results over all 
test conditions, the Press method seemed to yield the most accurate predictions* The 
largest error in most cases occurred at transonic speeds where shock-induced separa- 
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tion was encountered. In view of the difficulty in making accurate separation loca- 
tion predictions# RAXJET allows the separation location to be specified by the user 
as an option. 

Discriminatlgq streamline for separation region .- Ihe discriminating streamline 
concept is used to separate the reverse flo^^ region from the outer boundary-layer 
flow. The streamline shape is calculated by the discrete control volume analysis 
developed by Presz et al. (ref. 11). T.\e method solves integral forms of the con- 
tinuity and streanwise momentum equations to account for streamwise pressure gradi- 
ents and surface skin friction. An error- fonct ion- type velocity profile modified to 
allow reverse velocities is assumed. The solution is marched downstream from the 
separation point to the end of the afterbody. 

Discriminating streamline for jet wake region .- Two methods are provided for 
calculating the discriminating streamline shape in the jet wake region shown in fig- 
ure 2. The first is an integral method that accounts for entrainment effects based 
on the turbulent mixing analysis of Peters and Phares (ref. 20). The integral method 
assumes a three-parameter polynomial form for the velocity profile which is matched 
to the error-function profile of Press at the end of the afterbody. Integral forms 
of the continuity# axial-momentum, and turbulent kinetic energy equations are solved 
by a marching technique from the end of the afterbody to the reattachment point (at 
which the reverse velocity vanishes). The turbulent kinetic energy profile is 
related to the velocity profile by a mixing-length model for the eddy viscosity as in 
reference 13. 

The second method for treating the jet wake region is a simple extrapolation 
model. In this model# the slope of the discriminating streamline at the end of the 
afterbody (as calculated by the Press method) is linearly extrapolated to intersect 
the inviscid plume boundary. While this model is ad hoc# its results do approximate 
the results of the original conical discriminating streamline model used in earlier 
work by Presz (ref. 8). 


Mixing-Layer Solution 

The displacement-thickness distribution arising from entrainment into the jet 
mixing layer is calculated by the overlaid mixing model of reference 13 and computer 
program BOATAC (ref. 17). 60ATAC solves the parabolic mixing equations by a finite- 
difference marching procedure in cylindrical streamline coordinates. Turbulence is 
modeled by either a mixing-length or a two-equation (ke2) transport model. A laminar 
viscosity option is also included. The calculation is initialised at the nossle exit 
for attached flows or at the reattachment location for separated flows. A fixed 
number of grid points# evenly spaced in the stream function coordinate# span the 
mixing layer. This computational grid is overlaid onto the inviscid external and jet 
flow fields (also mapped into streamline coordinates) from which edge conditions and 
pressure gradients are determined. The solution is marched downstream a user- 
specified distance that is sufficient to account for the near-field interaction 
(typically 10 to 15 exit radii). An effective plume boundary which accounts for mass 
entrainment into the mixing layer and inviscid plume blockage is calculated as 
described in reference 14. Then the inviscid plume boundary is subtracted from this 
effective boundary to define a displacement thickness related only to entrainment. 
This displacement-thickness distribution is treated like an extension of the 
boundary-layer displacement distribution and is added to the inviscid plume boundary 
by underrelaxation in the iteration procedure described in the next section. 
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Viscous-Inviscid Iteration Procedure 


The component solutions are combined into the following iterative procedure 
(presented in fig. 3) : 

1* Calculate the inviscid external flow field over the effective body of 
revolution. 

2. Calculate the inviscid jet exhaust boundary and flow field. If the flow is 

separated# the inviscid jet exhaust solution is frozen after 10 iterations. 

3. Calculate the boundary-layer displacement thickness. 

4. Calculate the separation location if required. The predicted separation 

location may be frozen after a user-specified number of iterations. 

5. Calculate the shape of the discriminating streamline if required. Note that 

this calculation is repeated for every iteration if separation occurs# even 
if the separation location is frozen. 

6. Calculate the mixing-layer displacement thickness. 

7. Correct the effective body geometry by adding the displacement thickness and 

discriminating streamline to the original body and inviscid plume bound- 
ary. The relaxation procedure given in reference 14 is used to add 
these corrections. 

8. Repeat steps (1) to (7) for the desired number of iterations. 

At the end of step (8)# the user must manually check the solution for convergence 
(usually by monitoring surface pressure and displacement distributions). Note that 
the order of calculation above is different from that in reference 14# although a 
single iterative loop is still used. 


COMPARISONS BETWEEN PREDICTIONS AND EXPERIMENT 

Predictions of program RAXJET are compared in figures 4 to 6 with experimental 
afterbody pressures from references 21# 22# and 23. All predictions were made with 
77 by 39 grid points in the inviscid external region# 41 grid points in the mixing 
layer# and 41 grid points in the inviscid jet exhauf^t. The two-equation turbulence 
model was used for all mixing-layer calculations. For the separated flow cases# the 
experimental separation location and extrapolated jet wake model were used. 

The comparisons for attached boundary-layer flow on a circular-arc boattail 
nozzle with \/D * 1.768 and D /D ■ 0.51 are shown in figure 4 at several nozzle 
pressure ratios (NPR) and free-stream Mach numbers. Predicted and experimental pres- 
sure distributions agree very well in almost all cases. The only significant differ- 
ences occur at » 0.96 and NPK - 2.0 (fig. 4(c)) where shock— boundary-layer 
interactions are expected. 

The comparisons for slightly separated flow on a circular-arc nozzle with 
\/D « 1.0 and for significantly separated flow on a nozzle with \A> = 0.8 are 
given in figures 5 «nd 6# respectively. Again the agreement between prediction and 
experiment is excellent except for the i/D » 0.8 nozzle at Mach numbers of u.9 and 
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0.94 (fig. 6(c)). While converged solutions were obtained and reasonable levels for 
the separation plateau pressures were predicted, the results at Mach numbers of 0.9 
and 0.94 are quite poor in the region of shocks and in the jet wake. Furthermore, 
calculations at these Mach numbers made with the Presz model gave poor predictions of 
the separation location, and attempts to improve the predictions by iterating on the 
separation location sometimes led to divergent solutions. Thus it appears that 
strong interaction modeling is needed to treat flows with shock-induced separation. 
Therefore, RAXJET is not recommended when separation Mach numbers are supersonic. 

Calculations made with the Presz model to predict separation location are com- 
pared to those made with the experimentally determined location in figure 7 for the 
\/D - 1.0 and \/D = 0.8 nozzles. In both cases the integral wake model was used, 
and the predicted separation location was frozen after the first iteration. For the 
\/D =1.0 nozzle/ the separation was predicted at Ax/D = 0.71, well upstream of the 
experimental value of 0.84 which resulted in a significant underprediction of the 
pressures near the end of the afterbody. For the \/D = 0.8 nozzle the predicted 
separation was at Ax/D = 0.54, compared to an experimental value of 0.51, and much 
better agreement is seen. In principle, better predictions of separation location 
could be obtained by updating the predicted location for additional iterations since 
the inviscid pressure should be approaching a more correct solution. However, the 
Presz model was found to be so sensitive to the imposed pressure distribution that 
significant amounts of underrelaxation were required to achieve convergence, and the 
predicted separation locations were not found to be consistently improved over those 
obtained on the first iteration. This sensitivity of separation predictions is typi- 
cal of separation criteria that are based on calculating a pressure rise to separa- 
tion (see ref. 10). 

Predicted results with the integral and extrapolated jet wake models are com- 
pared in figure 8. Both sets of calculations used the experimental separation loca- 
tion. Slight differences are noted between the results of the two methods, and the 
judgment as to which is the better method seems to depend on which aspect of the 
prediction is most important to a particular problem. From these and similar results 
other flow conditions, the integral method seems to yield more accurate flow-field 
predictions, while the extrapolated method seems to yield more accurate drag 
predictions. 

Predicted and experimental boatti il pressure drag for the data of figures 4 
and 6 are compared in figure 9. Predicted drag obtained with the Presz separation 
model and the integral wake model and with the experimental separation location and 
the extrapolated wake model are shown for the separated flow cases. The drag calcu- 
lated with the experimental separation location and the integral wake model fell 
between these two curves. Predicted results from reference 6 which neglect entrain- 
ment effects downstream of reattachment are also included. Generally, the agreement 
between the present results and experiment is good with noticeable improvement over 
the results of reference 6. This improvement is primarily due to the inclusion of 
jet entrainment effects. Although there is an undesirable sensitivity of the present 
predictions to the separation and jet wake model used, the present method does cor- 
rectly predict the qualitative behavior for a given model. Thus, for parametric drag 
studies it is recommended that a particular combination of separation and wake models 
be used consistently. 

In figure 10, predicted pitot pressures across the jet and mixing layer are 
compared with the data of Mason and Putnam (ref. 24). Figure 10(a) gives the results 
for a nozzle with attached flow at a Mach number of 0.4 and a nearly fully expanded 
jet (NPR • 2.0). The agreement is reasonably good, except that the spread rate of 
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the mixing layer is underpredicted. Cline and Wilmoth (ref. 25), who used a Navier- 
Stokes solution method with essentially the same turbulence model, found a simi?.ar 
trend. The results of reference 14 suggest that at least some of the discrepancy may 
be due to deficiencies in the turbulence modeling. However, the overall afte'^h?» 
flow-field predictions were found to be relatively insensitive to the turbut*i .c~. 
model used in the mixing- layer calculation. 

Figure 10(b) gives results for a nozzle with separated flow at a Mach number of 
0.8 and a moderately underexpanded jet (NPR » 5.0). The ^^irediction was made with the 
integral jet wake model. Mote that the wake defect predicted 1:^ this model agreas 
with the experimental profile nearest the nozzle exit (Ax/r^ " 0.1). Also, the pre- 
dicted structure of the invlacid plume flow field is in excellent agreement with the 
experimental structure. The locaticn and strength of the barrel shock and Mach disk 
normal shock and the size of the Mach disk streamtube are all correctly predicted. 
Since this shock structure is known *^0 be sensitive to the external pressure distri- 
bution imposed (see ref. 24), the present results suggest that the basic plume inter- 
action model is reasonably correct. 


DESCRIPTION OP COMPUTER PROGRAM 

The organization of computer program RAX JET is presented in figure 11. The 
program is written in overlay form and consists of the main overlay (0,0) and six 
primary overlays (fig. 11 (a)). The main overlay, RAXJET, controls the input and 
output of data and the viscous-inviscid iteration scheme. Primary overlay (1,0) 
contains the program RAXBOD which calculates the inviscid external flow. Primary 
overlay (2,0) contains the progreun SCIPAC which calculates the inviscid jet exhaust 
flow. Primary ov''rlay (3,0) contains the program VISCOUS which calculates the 
boundary-layer growth and performs the separation analysis. Primary overlay (4,0) 
contains the program BOATAC which performs the overlaid mixing analysis. Primary 
overlay (5,0) contains the progr^un INPT which reads the input data, and prlioary 
overlay (6,0) contains the program OUTPT which prints the computed results. 

The program uses six disk files. Input data are obtained from file TAPES and 
the results are written to file TAPE6 which is equated to file OUTPUT. The remaining 
disk files are used internally by the program and are described later. RAXJET 
requires about 76 000 octal storage locations on the Control Data CYBER 175 con^uter 
and executes 20 viscous-inviscid iterations in approximately 12 minutes of CPU time 
for separated flow on an afterbody with a fully expanded jet. Attached flow compu- 
tations generally require fewer iterations and therefore less conqputer time. Under- 
expanded or overexpanded jets generally require more computer time. 

A brief description of each of the main routines is given in the following listi 

1. Program RAXJET .- RAXJET initiates a run by calling program INPT to read the 
input data. Then an initial call is made to program RAXBOD to set up the inviscid 
external grid. This grid is used to define various crid indices used by other rou- 
tines and to Initialize various geometric arrays. RAXJET then starts the viscous- 
inviscid loop Iv calling RAXBOD, SCIPAC, VISCOUS, and BOATAC, successively. Then 
program OUTPT is called to print the results. Finally, RAXJET calculates the new 
effective body and repeats the viscous-inviscid iteration until the specified nun^r 
of iterations is consisted. RAXJET also calls subroutine AINTPL to Interpolate the 
original body geometry onto '>.e Inviscid grid (which changes sll^tly during each 
iteration) and subroutine SMOOTH to smooth the effective body geometry for separated 
flows. 
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2. Program RAXBOD >" RAXBOO solves the exact nonlinear potential flow equation 
for the inviscid external flow and consists of the primary overlay (1/0) and four 
secondary overlays (see fig. 11(b)). RAXBOD first calls ONE1 (overlay (1/1)) to set 
up the tangential grid and ONE2 (overlay (1,2)) to set up the normal grid. Since the 
effective body shape may change with each iteration, these routines must be called 
before each inviscid calculation. For the first iteration only, subroutine ESTIM is 
called to initialize the perturbation potential to zero. For each succeeding itera- 
tion, potential values from the previous iteration are used as the initial estimate. 
Program ONE3 (overlay (1,3)) is then called to solve the exact potential equation by 
the South- Jameson relaxation scheme. The relaxation is continued until either the 
specified maximum number of relaxation iterations are perforu^ed or the specified 
convergence tolerance is reached. Finally, program ONE4 (overlay (1,4)) is called to 
calculate the flow quantities (velocity, pressure, etc.) from the perturbacion poten- 
tial and to transfer the results to disk file TAPE4. 

3. Program SCIPAC .- SCIPAC solves the inviscid jet exhaust region by a shock- 
capturing — shock-fitting procedure and consists of the main procjrain, 24 subroutines, 
and 6 function i-ubroutines (see fig. 11(c)). SCIPAC initializes a jet calculation by 
calling SCIPPB v,hich reads the external flow properties from TAPE4 and determines the 
pressure to be specified along the plume boundary. Then SCIPIN is called to generate 
the initial conditions at the nozzle exit from the input data and to discretize the 
initial flow conditions for the finite-difference marching procedure. Also, from 
nozzle exit conditions, SCIPIN estimates the length of the first inviscid cell *A’bich 
is used to terminate the marching procedure. Then SCIPCT is called to integrate the 
finite-difference flow equations and to perform the iterative Mach disk calculation. 
For the Mach disk calculation, disk files TAPEl and TAPE2 store conditions upstream 
of the predicted Mach disk location while trial solutions downstreaiA are attempted. 
Details of the Mach disk procedure are given in reference 17. Depending on certain 
input param3ters and the predicted inviscid cell length, the marching procedure 
terminates either downstream of the Mach disk or at the end of the first or second 
inviscid cell. Properties downstream of this termination point are calculated b\' a 
simple isentropic expansion to local external pressure by calling subroutine SCIPDK. 
Jet exhaust properties at axial locations corresponding to the external inviscid grid 
are transferred to disk file TAPE3 by calling SCIPOT. Finally, SCtPPL is called to 
determine the coordinates of the inviscid plume boundary and to store them for later 
use. 


^^ ^ gam VISCOUS .- Program VISCOUS calculates the boundary-layer growth, 
separation location, and discriminating streamline shape calls to subroutines ONE, 
TWO, and THREE (tiee fig. 11(d)). Subroutine ONE calculates the boundary-layer dis- 
placement thickness, skin friction, and jower-law exponent by a modif.^ed Reshotko- 
Tucker integral procedure with the invijcid pressure gradient imposed by the external 
flow solution. Subroutine TWO calculates the separation location by the Press con- 
trol volume method. Subroutine THREE calls B834 which calculates the difi criminating 
streamline shape by the modified discrete control volume method of Press et al. If 
the Integral wake model is specified, B834 calls JET to calcvlate the discriminating 
streamline from the end of the body to the end of the recirculation region. If tne 
integral wake model is not specified, B834 linearly extrapolates the discriminating 
streamline until it intersects the inviscid plume boundary. 

5. Pr ogram BOATAC .- Program BOATAC calculates the turbulent mixing layer and the 
plume displacement thickness and consists of the primary overlay (4,0) and two 
secondary o^'erlays (see fig. 11(e)). BOATAC first calls BOATSA (overlay (4,1)) which 
in turn calls (1) BOATIF to read the inviscid flow field ivorx TAPE3 and TAPE4, 

(2) BOATIP and BOATKE to generate the initial velocity, tenq>erature, and turbulence 
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profiles/ and (3) BOATSB to map the profiles into equally spaced streamline coordi- 
nates. BOATAC then calls BOATMl (overlay (4/2)) which in turns calls (1) B0ATS3 to 
integrate the parabolized flow equations/ (2) BOATDS to calculate the plume dis- 
placement thickness, (3) BOATEN to estimate the necessary computational boundary 
growth rate, (*") B0ATM2 to reset the discretized profile to btgin the next marching 
step/ and (5) BOATOT to transfer the results to disk file TAPE2 at axial locations 
corresponding to the external grid. BOATAC starts at the nozzle exit for attached 
flows or at the reattachment location for separated flows and marches downstream the 
user-specified number of jet radii. Mixing-layer properties downstream of this axial 
location are not calculated. 

6. Program INPl INPT reads the innut data from TAPES, checks for possible 
input errors, and prints the appropriate warning messages (fig. 11(f)). If certain 
RAXBOD grid parameters are not input, INPT calls subroutine OPTIM which attempts to 
optimize those parameters not given so that a maximum number of grid points are 
located on the afterbody consistent with an accurate solution of the total flow 
field. Finally, INPT calls subroutine lOUTPT to print the input data. 

Program OUTPT .- OUTPT prints the surface quantities, drag coefficients, and 
convergence history (fig. 11(f)). If requested, OUTPT calls subroutine OPFOUT to 
generate and print a composite viscous-inviscid flow field. 


Description of Input Data 

The program input data are entered by a title card followed by eight namelists. 
A sample input data deck is given in the appendix. Some of the geometric input 
quantities are defined in figure 12, and certain input parameters for the off-body 
flow- field output are given in figure 13. 

1. Title card .- The first card nf each data deck consists of 80 alphanumeric 
characters that identify the case* The eight namelists must then appear in the 
following order. 

2. Namelist CNTRLN .- CNTRLN contains the input parameters that control the main 
logic of the program: 

INAX Integer that specifies the number of grid points in the streamwise 

direction for the inviscid extarnal flow solution (81 maxim im). 

The default value is 81. 

JMAX Integer that specifies the number of grid points in the body normal or 

radial direction for the inviscid external flow solution (41 maximum). 
The default value Js 41. 

MIT Integer that specifies the r*imber of relaxation cycles used per viscous- 

inviscid iteration for the x.iviscid external flow solution. The 
default value is 20. 

MVI Integer that specifies the number of viscous-inviscid iterations 

(30 maximum). The default value is 1. 

ISBP Integer that specifies the separation model used: 

If ISEP « 0, no separation model is used. 

If ISEP ■ 1, the Press control volume method is used. 
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If ISEP = 2, the separation location is specified by the user. 

The default value is 0. 

NSEP Integer that specifies the number of viscous-inviscid iterations after 

which the separation location is held fixed (needed only if ISEP = 1), 
The default value is 1. 

IWAKE Integer that specifies the jet wake model (needed for I.^EP - 1 or 2): 

If IWAKE « 0, the integral wake i^odel is used. 

If IWAKE ==1/ the extrapolated wake model is used. 

The default value is 0. 

XSEP Axial coordinate of the separation location nondiraensionalized by REEL 

(needed only if ISEP « 2). The default value is XM. 

lUI Integer that specifies the type of units used for input quantities: 

If lUI = 0, lengths are input in feet/ temperatures in Rankine/ and 
pressures in pounds per square foot (psf). 

If lUI = 1/ lengths are input in meters, temperatures in kelviU/ and 
pressures in atmospheres (1 atm = 101.325 kP?). 

The default value is 0. 

lUO Same as lUI except it applies to output quantities. 

3. Namelist ESC .- ESC contains the free-stream conditions and reference 
quantities: 

GAM Ratio of specific heats. The default value is 1.4. 

AMINE Eree-stream Mach number. No default value is specified. 

PT Eree-stream total pressure in pounds per square foot or atmospheres. 

The default value is 2116 psf or 1 atm, depending on lUI. 

TT Eree-stream total temperature in ®R or K# The default value is 530® R 

or 294 K, depending on lUI. 

REEL Reference length used to convert input coordinates to feet (lUI » 0) or 

meters (lUI = 1). All other inputs having dimensions of lengths are 
assumed to be nondimensionalized by REEL. The default value is 1.0. 

SREF Reference area in square feet or square meters used to calculate drag 

coefficients. The default value is (ir/4) (REEL)^. 

ALPHAS Array contaiiiing the mole fractions of the six species in the following 
order: N 2 / O 2 / CO 2 / H 2 O/ CO, and CH^. The default values are the 

nominal values for air: ALPHAB(I) « 0.79, ALPHAE(2) » 0.21, 

ATiPHAE(3) » ALPHAE(4) « ALPHAE(5) » ALPHAE(6) - 0.0 

4. Namelist JETDAT .- JETDAT contains the* nozzle exit conditions and control 
parameters for the inviscid jet exhaust solution: 

XNPH Nozzle pressure ratio, ratio of jet total pressure to free-stream static 

pressure. No default value is specified. 
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EMJET Nozzle exit Mach number (EMJEP > 1.05). The default value is 1.05. 

TTJET Nozzle exit total temperature in ®R or K. The default value is 530® R 

or 294 K, depending on lUI. 

RJET Nozzle exit radius nondimensionalized by REEL. No default value is 

specified. (See fig. 12.) 

THLIP Nozzle exit lip angle in degrees. The default value is 0.0. (See 
fig. 12.) 

NMAXJ Integer that specifies the number of radial grid points used in the 
inviscid jet solution {81 maximum). The default value is 81. 

IGAS Integer that specifies the thermodynamic option for the inviscid jet 

solution: 

If IGAS 0, a perfect gas is assumed and GAMJ » 1.4. 

If IGAS = 1^ a perfect gas is assumed and GAHI must be input. 

If IGAS =* 2, a calorically imperfect gas is assumed and GAMJ will be 

calculated. 

If IGAS = 3, a perfect gas is assumed but GAMJ will be calculated 
based on nozzle exit conditions and held constant downstream. 

The default value is 0. 

GAMJ Ratio of specific heats for inviscid jet exhaust. See IGAS for defaults. 

lOK Integer that controls the solution method downstream of a Mach disk: 

If IDK « 0, an isentropic decay Is imposed. 

If IDK 1, a finite-difference shock-capturing solution will be 
attempted. 

The default value is 0. 

ALPHAJ Array containing the mole fractions of the six jet species. Same as 
ALPHAE (namelist FSC) except it applies to the inviscid jet. 

5. Namelist MIXDAT .- MIXDAT contains the input quantities that control the over- 
laid mixing-layer calculation: 

IVIS Integer that specifies the turbulence model. 

If IVIS « -2, the ke2 model is used. 

If IVIS * 0, the Prandtl mixing-length model is used. 

If IVIS • 2/ a laminar model is used. 

The default value is -2. 

MPSI Integer that specifies the number of radial grid points used in the 

mixing-layer calculation (50 maximum). The default value is 31. 

FDL Factor that multiplies the marching step size (FOL < 1.0). The default 

value is 1.0. 

FFF Mixing-length scale used in the Prandtl mixing model in the near-field 

region. The default value is 0.065. 

GOG Mixing-length scale used in the Prandtl mixing model in the fully 

developed region. The default value is 0.080. 


It 
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SIGMA Tiirbul«nt Prandtl number. The default value is 1.0. 

TCONT Maximum allowable temperature change in one marching step in ®R or 

depending on lUI. The default value is 400/(MPSI - 1) K. 

XJENT Axial distance in jet radii from the nozzle exit where the mixing 
calculation will be terminated. The default value la 15. 

KMAXJ Integer that specifies the number of mapped radial grid points in the 

inviscid jet solution that are used In the overlaid mixing calculation 
(25 maximum). The default value is 25. 

KMAXE Same as KMAXJ except it applies to the inviscid external flow solution. 

TKEJ Factor that multiplies the initial turbulent kinetic energy in the jet 

at the nozzle exit. The default value is 1.0. 

TKEX Same as TKEJ except it applies to the external flow side of the initial 

profile. 

I^Amelist REIiAX .- RELAX contains various relaxation parameters used in the 
inviscid external flow solution and vlscous-invlscid Interaction: 


RP1 Relaxation factor for subsonic points in the inviscid external flow 

solution (0 < RF1 < 2). The default value is 1.4. 

QF3 Supersonic damping factor for Improving Iterative stability in the 

Inviscid external flow solution (QF3 > 0). The default value is 0.1. 


COVERG Convergence criterion for inviscid external flow solution. Inviscid 

relaxation cycles stop for a given iteration when the average residual 
Is less than COVERG/(IMAX - 1)^ or MIT cycles are completed, whichever 
occurs first. The default value is 1.0. 


Wn and Underrelaxation parameters for adding the displacement thickness and 
WBP discriminating streamline to the body. WB is the value for the first 

iteration and WBF Is the value for the final iteration. For inter- 
mediate Iterations, the underrelaxation factor Is varied linearly 
from WB to WBF. (WB and WBF must be in the range from 0 to 1.) The 
default values are 0.5. Larger values of WB or WBF give more 
under relaxation. 


^Amelist GRID .- GRID contains the parameters that control the body-fitted 
grid used in the inviscid external flow solution: 

XBT Axial location of the start of the afterbody nondimenslonallzed by REFL. 

No default value is specified. (See fig. 12.) 

XM Axial location of the end of the body nondimensionalized by REFL. No 

default value is specified. (See fig. 12.) 

DXIDXO Derivative of the tangential coordinate stretching function at the nose, 
(drydX)^.Q. A reasonable value can be estimated from 
DXIDXO • (IMAX - DAIC where Af, is the desired grid size at the 
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nosot If no inpvit value is given# the program will attempt to 
calculate an optimum value. 

DXIDXM Derivative of the tangential stretching fvmction at the end of the body, 
(dr/dX)jj.jjl^|. A reasonable value can be estimated from 
DXIDXN •• (IMAX - ^)Ar^ where is the desired grid size at the 

end of the body. If no input value is given, the program will attempt 
to calculate an optimum value. 

CXM Fraction of streamwise grid points that lie on the body (0 < X < XM). 

(CXM must be in the range from 0 to 1.) The default value is 0#75. 

Note that if no input value of DXIDXM is given, CXM may be changed in 
an attempt to find an optimum value for DXIDXM. 

DNDYO Derivative of normal coordinate stretching function at the body svirface, 
(dn/dY) . A reasonable value may be estimated from 
DNDYO AV)^ An /AY where An is the desired grid size at the 

surface, a is the stretching exponent, and AY » 1/(JMAX - 1). If no 
input value is given, the program will attempt to calculate an optimum 
value. 

ALF Exponent in the normal coordinate stretching function. (ALF > 0.) The 

default value is 1.3. Larger values of ALF move the last finite 
value of n farther away from the body and smaller values move it 
closer. 

Namelist OUTPTC .- OUTPTC contains parameters which control the printing of 
results. Some of the parameters used to generate the composite flow-field output are 
defined graphically in figure 13. Figure 14 shows how the composite flow field is 
constructed. 

TOUT Integer that specifies the iteration number at which the printing of 

results will begin. All Iterations from lOUT to MVl will be printed. 
The default value is 0. 

I AFT Integer that controls the amount of results printed: 

If lAFT >■ 0, results will be printed starting at the nose* 

If lAFT - 1, results will be printed starting at the first grid point 

ahead of the afterbody. 

The default value is 0. 

lOFF Integer that controls the printing of the off-body flow-field results: 

If lOFF " 0, no off-body results will be printed. 

If lOFF - 1, off-bod>/ results will be printed. 

The default value is i). 

XSTRT Axial location (nondim<inaionalissd by RfiFL) at which the printing of off- 
body results will begin* No default value is specified. 

XEND Axial location (nondimenaionalised by RBFL) at which the printing of off- 

body results will end. No default value is specified. 

ROB Radial outer boundary (nondimenaionallsed by RSFL) for printing off-body 

results. No default value le specified. 
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IMAP Integer that specifie;i the number of mapped output points between the 

body (or inviscid p .ume boundary) and r “ ROB (see fig. 13). 

(IMAP < 100 - JMAP.) The default value is 25. 

JMAP Integer that specifies the number of mapped output points across the jet 

exhaust (0 < r ^ Inviscid plume boundary). (JMAP < 100 - IMAP.) 

The default value is 25. 

ALFO Exponent in radial stretching function used to distribute output points 

in off-body results. AFLO > 1 will pack points in the boundary layer 
and mixing layer. ALFO * 1 will give equally spaced points. The 
default value is 1. 


9. Namelist GEOMN .- GEOMN contains the body geometry data. Note that a finite- 
length forebody geometry must be input even if only afterbody results are desired. 
The forebody nose may be blunt or pointed, and the input coordinates are assumed to 
be referenced to the nose (x « 0, r * 0 in fig. 12). 


IXY Integer that specifies the number of body input coordinates 

(140 maximum). No default value is specified. 

lORDER Integer that specifies the type of interpolation used to determine 
the geometry at the inviscid grid points: 

If tORDER * 1, first order will be used. 

If lOROER » 2, second order will be used. 

The default value is 1. 


XO Array of IXY axial body coordinates (nondimensionalized by REPL) with 

origin at the nose of the body of revolution. No default values are 
specified. 

YO Array of IXY radial body coordinates (nondimensionallzed by REFL) with 

origin at the nose of the body of revolution. No default values 
are specified. Note that the first point of the XO and YO arrays must 
correspond to the nose, XO(1) » 0, and che last point to the end of 
the body, XO(IXY) * XM. 


Description of Output 

An example of the printed output is given in the appendix for the sample input 
data a4.so presented there. The first two pages are duplicate output with RAXJET 
bannerised followed by the case title, date, and time of run. Warning messages gen- 
erated by the input data are printed on the third page and diagnostic messages from 
the grid optimization routine are printed on the fourth page. The fifth page con- 
tains the program title and abstract, case title, and free-stream and jet exhaust 
conditions. On the sixth page, the gas composition and remaining input parameters 
are printed. The next two pages list the input geometry. The tabular Input geometry 
has beer, augmented with points that define a cylindrical sting starting at the end 
of th ^ body as an initial guess for the plume shape. 


The following pages list the results. For each iteration starting with lOUT the 
case title, Iteration number, free-stream Mach number MO, jet pressure ratio, and 
reference length L are given. If boundary-layer separation is predicted, the 
separation and reattachment locations are also printed. A tabulation is given of 
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Axial grid point number I 
Body axial coordinate X/L 
Body radial coordinate H/b 
Pressure coefficient CP 
Local Mach number ML 
Skin-friction coefficient CF 

Body radius corrected for the discriminating streamline and inviscid plume 
shape RDS/L 

Effective body radiuc REFF/L 

Radial location of the outer edge of the viscous layer RVI/L 
Displacement thickness DEL*/L 
Boundary-layer momentum thickness THETA/L 
Compressible shape factor H 

Note that REFF/L is the effective body obtained by underrelaxation, while 
RDS/L + DEL*/L is the calculated displacement body for a given iteration. Note also 
that the definition of DEL^/L for the plume mixing layer used here differs from that 
presented in reference 14. Here, values of DEL*/L for X/L downstream of the nozzle 
exit include only the contribution due to jet enrrainment and do not include the 
inviscid plume blockage. However, the effective body definition used here is the 
same as in reference 14. After the final iteration is printed, the following drag 
coefficients for each iteration are tabulated: 

Afterbody pressure drag CDP,AFT 

Afterbody friction drag CDT,AFT 

Total (pressure and friction) afterbody drag CDT,AFT 

External pressure drag for the complete body of revolution (forebody and after- 
body) CDP,BOD 

Body friction drag CDF, BOD 

Total body drag CDT,BOD 

Next# the maximum residual RPNAX and average residual RPAV6 for the RAXBOD calcula- 
tion are listed as an indication of the convergence history. 

If off-body flow-field results are requested# the remaining pages tabulate 

Axial location X/L 

Radial location R/L 
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Nondimens lonalized pressure P/PINF 
Nondimens lonalized temperature T/TINF 
Nondimens ionalized anial velocity U/UINF 
Nondimensionalized radial velocity V/UINF 
Mach number M 

Ratio of specific heats GAMMA 

Nondimensionalized pitot pressure PT/PTINF 

Pressures, temperatures/ and velocities are determined by linear interpolation 
directly from the appropriate flow solution as indicated in figure 14 except for 
points in the boundary layer (denoted by BL) . For the BL points, edge conditions and 
pressures are determined by interpolation from the inviscid solution, streamwise 
velocities are calculated from a power-law profile, and temperatures are calculated 
by assuming a constant total temperature across the layer* Axial and radial velocity 
components are determined from the streamwise values by assuming a linear variation 
in streamline slope across the layer from the local body slope at the surface to the 
inviscid streamline slope at the outer BL edge. For all points, local Mach number 
and pitot pressure are computed from these interpolated and calculated quantities* 
Note that for supersonic points, the pitot pressure is obtained by normal shock rela- 
tions for an adiabatic, perfect gas. Points in the inviscid jet region that are 
computed from isentropic theory are denoted by (*) and points in the mixing layer 
downstream of where the mixing solution terminates are denoted by (♦♦)• Note that 
the axial location of points in the inviscid jet and mixing layer may not correspond 
exactly to the grid location in the external flow* These points correspond to actual 
marching step locations that are nearest the external grid location* No attempt was 
made to interpolate these points in the axial direction since the discrepancy is 
usually (juite small. Also note that flow-field results are not given within the 
separated region. 


Instructions to the User 

The following special instructions are provided to assist the user in applica- 
tions of program RAXJET. These instructions do not cover all possible user problems 
but are intended to assist the user in anticipating cases where the physical and 
computational models may not give accurate results. Many of the problems associated 
with inappropriate input data are flagged by the program followed by explanatory 
warning messages, and these are not discussed here. 

Grid generation .- The user has input control over the streamwise and normal grid 
used by RAXBOD, and the radial grid used by SCXPAC and BOATAC. The streamwise grid 
used by VISCOUS for the boundary- layer calculation is the same as that for RAXBOD* 

For most problems, the built-in procedures for optimizing the RAXBOD grid parameters 
give satisfactory results as long as the afterbody length comprises approximately 
10 to 25 percent of the total body length* For shorter afterbodies, other unusual 
geometries, or much fewer than the maximum allowable number of grid points 
(XMAX,JMAX), the user should verify that the grid is adequate. For the SCIPAC calcu- 
lation, accurate shock-capturing may require nearly the maximum allowable number of 
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points (NMAXJ ==81)* If fewer points are used, SCIPAC gives less accurately defined 
shocks and may encounter problems when computing cases with small Mach disks or regu- 
lar barrel shock reflections or when attempting to extend the computation beyond the 
first Mach disk or reflection (IDK =1). If such problems are encountered, the user 
should try increasing NMAXJ. For the BOATAC calculation, the default value of 
MPSI * 31 may be adequate for the basic entrainment interaction; however, larger 
values of MPSI are recommended if accurate mixing-layer profiles are desired. 

Convergence? .- The convergence history of the RAXBOD calculation is given as part 
of the output. If convergence is not achieved 1^ the final iteration, the number of 
global iterations MVI, the number of relaxation cycles MIT, or both should be 
increased. An indication of the global convergence of the viscous-inviscid calcula- 
tion can be obtained by monitoring the afterbody pressure drag (CDP,AFT) or selected 
pressures and, for separated flows, the separation location (if ISEP =1, 

NSEP > 1). 

Attached flows (ISEP = 0) .- Calculations in which attached flow is anticipated 
generally converge in 5 to 10 iterations with most input parameters set at their 
default values. Plows with incipient separation may also be ru.* as attached flows 
but may require greater amounts of underrelaxation {increased wB and WBF) and more 
iterations to converge. 

Separated flows (ISEP > 0) .- Separated flows require considerably more user 
intervention in selecting appropriate control parameters than attached flows. The 
greatest source of error generally lies in prediction of the separation location. 

More accur«'\te flow predictions can be made if the separation location is known and is 
input by the user (ISEP =2). If the separation location is not known, it is recom- 
mended that a solution first be attempted as for attached flow (ISEP « 0). If the 
solution diverges or if the boundary-layer shape factor exceeds a value of 2, then 
assume that the flow is separated and attempt a solution with ISEP * 1. A value 
of NSEP = 1 (the default value) is recommended which bases the separation location 
prediction on the initial inviscid solution. This procedure yields reasonable 
results except where the flow is approaching shock-induced separation. RAX JET is not 
recommended for strong shock-induced separation. More accurate predictions of the 
separation location can sometimes be obtained by allowing the predicted location to 
be updated over several iterations (NSEP > 1), but larger amounts of underrelaxation 
are usually necessary (WB = 0.75, 0.5 < WBF < 0.75, typically). The cl oice of the 

wake model depends on the type of result that is of primary interest. If aft^>“body 
pressures and drag are the main interest, IWAKE » 1 is recommended. If mixing- 
layer profile results are important, IWAKE « 0 is recommended. 

Real gas thermodynamics .- The real gas thermodynamics used in SCIPAC and BOATAC 
contain built-in data for the six species, N 2 # O 2 / CO 2 / H 2 O, CO, and CH^. Alternate 
species may be used by replacing the appropriate polynomial coefficient data in sub- 
routines SCIPCF and BOATSD with that for the alternate species. Detailed definition 
of these coefficients is given in reference 17. 


CONCLUDING REMARKS 

A viscous-inviscid interacton method has been developed to compute the subsonic 
and transonic flow over nozzle afterbodies with supersonic jet exhausts. The method 
accounts for the interaction between the inviscid external flow, an underexpanded jet 
exhaust, an attached or separated boundary layer, and the turbulent jet mixing layer. 
A computer algorithm called RAXJBT has been written which iteratively combines 


19 



solutions for the various flow regions by a relaxation procedure* The results of the 
method are in good agreement with experiment for afterbodies with attached flow up to 
moderately high transonic speeds* For afterbodies with separated flow, the agreement 
is good except when shock-induced separation is encountered* The results indicate a 
need for improved modeling of separated flows where strong interactions occur. The 
method correctly predicts the effects of jet entrainment and plume blockage on nozzle 
drag. Good agreement is obtained with flow-field data for moderately underexpanded 
jet exhausts* 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 4, 1982 
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APPENDIX 


SAMPLE INPUT AND OUTPUT FOR PROGRAM RAXJET 


Sample Input Data 


U/n*y,tJ» ClHCUlArt A^iC HiuTTAU • s nPK s ^,0 

*CmTHLA' I«AAs7;,J'’Axsl9# ''Vlsi> 0 ,t-'lTrU 0 , ISfcHsUlUOsl > 

s.rsc A'"lfjFru.tt»PTs^l 10,, T rsb60,,«E.f L* 0 ,b * 

»Jt TUAT T TJ^ TsbiO, » WJt Tso,i?S,ivMAXJs9 1 » IGAS*2, KJKSOi 

fcMjeTst,()S % 

»<^IXIJAT MPSIsul b 

»M|-:LAX ■^fciau.TS, NHb so,S :b 

IGklU txrs,/S, *rtlSH, 0 ,XHsO,« s 

JU.jTHTC li»Ulsi»u,lAFTsl,lUFFsi,X5TWTsd,w,XtNUs9,a,ROHso,7S 4 
iLEOMiM 1XVS09, 

xO(l)*i>,U»#l,,2,#.^#,u,,b,,b»,7,»»8»,^»l*i)»l»l»l»2#l»4#l,tt,l *W4#1,7 


b,?#b, S»*5,A,S,S,^i,b»S,7,5,^»5,^fb,o,b, 1 ,b,2fb,3,b,4#b,Sfb,bfb,7 ^6,^? 
b, 9 , 7 , 0 , 7 .l, 7 . 2 , 7 , 3 , 7 ,U, 7 . 5 , 7 ,O, 7 . 7 . 7 ,H, 7 . 9 ,b, 0 ,b,l,», 2 , 8 . 3 , 6 ,tt,R,Si 


« 


7 #b,bf 

7(.'(nsw»0,,02Sui3,,osa(jPh,,i>7bu3b,,i('ou*>;»,12bue9,,lboo77#|l75090, 
,200lu2, ,22!>1 lb, ,?b«i 2a, ,?7biut , , ioo 1 5« , , 32b l oo, , 3b0 1 79, , 375 1 92, 
,bt)02iib,,«2b21«.#‘<bii23i/,,<47 3oii»,,AH97ou,,«9a3bo,59*,5,,«96<<9o, 
,«8b932, ,<40bl uo, ,bo9oao, , 30 789 7, ,316759, ,255 4 
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Figure 4.- P'^edicteti and experimental afterbody pressures for attached flow 

on a circular-arc boattail nozzle with \/D * 1*768 and » 0.51* 
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Figure 4*- Continued. 
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Figure 11.- Organizational chart for progreun RAXJET 
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(b) Program RAXBOO 
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